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TRANSIENT SOLIDIFICATION OF A FLOWING LIQUID ON A COLD PLATE 
INCLUDING HEAT CAPACITIES OF FROZEN LAYER AND PLATE* 
by Robert Siegel and Joseph M. Savino 
Lewis Research Center 

SUMMARY 

An analysis is made of solidification within a flowing medium where a convective 
boundary condition is imposed on the moving interface of the solidifying layer. The con- 
figuration is one in which a frozen layer forms in a warm liquid as it is chilled while 
flowing over a plane wall that is convectively cooled on the opposite side. The analysis 
includes the heat capacities of both the wall and frozen layer. The solution for frozen 
layer thickness variation with time is given as a closed-form algebraic equation and in 
graphical form so that results can be readily evaluated for use in practical applica- 
tions. Some illustrative examples are given for flowing water in contact with Inconel 
plates of various thicknesses to demonstrate the effect of the heat capacities and the liq- 
uid and coolant temperatures on the growth of the ice layer and the temperature distribu- 
tions in the ice and wall. 


INTRODUCTION 

This report is concerned with the analysis of solidification when a convective heat 
transfer condition is imposed at a moving interface of a frozen layer. This type of 
boundary condition is encountered in important applications such as the solidifying of 
metal castings in molds, freezing of rivers, and solidification within liquid flow heat ex- 
changers that use a cryogen as the coolant. As discussed in references 1 and 2, rela- 
tively little analytical work and even less experimental work has been done for conditions 
where solidification is occurring within a flowing liquid phase. 

The specific configuration considered here is the solidification of a warm liquid as 
it flows over one side of a plane wall that is convectively cooled on the opposite side. 
Initially there is no coolant flowing, and the wall and liquid are at a constant tempera- 
ture that is above the fusion temperature of the liquid. Then a coolant is introduced 
whose temperature is below the freezing temperature of the liquid. As the wall cools 

*This material was presented at the ASME Winter Annual Meeting in Pittsburg, Pa. 
Nov. 17, 1967 as Paper No. 67-WA/HT-34. 



below the fusion temperature, the flowing liquid in contact with the plate begins to solid- 
ify. 

During the transient growth period, the heat supplied to the liquid - frozen layer in- 
terface is made up of the convection from the warm liquid and the latent heat of fusion 
generated by the freezing process. At steady state only the convective heat flux is sup- 
plied to the interface. The thickness of the steady-state layer is determined by the heat 
balance that this convective flux is exactly equal to the heat conducted through the frozen 
layer and wall and convected to the coolant. 

In references 1 and 2, we investigated the transient solidification of a warm liquid 
on a thin wall. The heat capacity of the frozen layer was accounted for, but the wall heat 
capacity was considered to be negligible. The analytical solution was adequate for ap- 
plications where a thick frozen layer forms on a thin wall. For thick walls, however, 
the wall heat capacity may be large compared with that of the frozen layer and may sig- 
nificantly retard the frozen layer growth. It is the purpose herein to develop a method 
that can be easily used to predict frozen layer growth when both the heat capacities of 
the frozen layer and the wall are accounted for. 

To obtain the frozen layer growth, the transient heat conduction equation was inte- 
grated within the frozen layer and within the wall. This yielded coupled integral equa- 
tions for the temperature distributions in the two media, and an integral equation for the 
frozen layer growth; this set of equations was solved iteratively. The frozen layer 
growth with time was obtained as a closed-form analytical equation that is easily evalu- 
ated for practical applications. Part of the growth relation is presented in graphical 
form so that hand calculations can be readily made. 

To illustrate the influence of the wall and frozen layer heat capacities on the solidi- 
fication rate, several examples are given for freezing of water. The examples illustrate 
the effect of varying the wall thickness, the coolant temperature, and the water tempera- 
ture. The transient temperature distributions in the ice layer and wall are also shown 
for a few example cases. 


SYMBOLS 


a thickness of cooled wall 
a’ dimensionless wall thickness, a/X., 
c specific heat 

ir 

F frozen layer capacity term in growth equation, see eq. (26) 
G integral terms defined by eq. (17b) 
h convective heat-transfer coefficient 
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I integral terms defined by eq. (16c) 

I integral terms defined by eq. (16d) 

J integral terms defined by eq. (16e) 

J w integral terms defined by eq. (16f) 

k thermal conductivity of solidified material 
thermal conductivity of wall material 
L latent heat of fusion 

q heat flux in x direction 

R dimensionless parameter, (X g /k)/|(l/h c ) + (a/k^J 

S dimensionless parameter, c p (T f - T c )/L 

T temperature 

T* dimensionless temperature, (T - Tj)/(T C - Tj) 

T' w dimensionless temperature, (T - T f )/(T C - T f ) 

W wall capacity term in growth equation, see eq. (26) 

X thickness of frozen layer 

x' dimensionless thickness of frozen layer, X/X 

s 

X g thickness of frozen layer at steady state 

x position coordinate measured from frozen layer- wall interface 
x' dimensionless coordinate, x/X g 

a thermal diffusivity, k/pc 

/3 dimensionless coordinate, x/a 

y dimensionless parameter, 1 + k^/hpa 

£ dimensionless coordinate, x/X 

r time from start of solidification 

t time during wall temperature transient preceeding solidification 

t' dimensionless time, rh^(T^ - Tj)/pLX g 

$ function of y defined in eq. (27) 

Subscripts: 
c coolant 

f at freezing temperature 



I liquid phase of solidifying substance 

s steady state 

w wall 

1 wall boundary in contact with coolant 

2 wall boundary in contact with solidifying material 

I, n successive iterative approximations 

r 


ANALYSIS 

The model used for this study is the one- dimensional configuration shown in figure 1. 

The x coordinate has its origin at the interface of the frozen layer and wall, and the 
wall extends to x = -a. A warm liquid at a fixed temperature T^ flows over one side 
of the wall, and it is assumed that the convective heat-transfer coefficient h^ is a con- 
stant. A transient solidification process can then be initiated by introducing a flowing 
coolant on the other side of the wall. It is assumed that the coolant is at a fixed tem- 
perature T c and provides a constant heat-transfer coefficient h c . After introduction 
of the coolant, the wall cools until the freezing temperature is reached on the surface of 
the wall (x = 0) exposed to the warm liquid. At this instant (r = 0, X(r) = 0) solidifica- 
tion is assumed to begin. The liquid - frozen layer interface is assumed to always be at 
the equilibrium fusion temperature T^. 

During the solidification process latent heat of fusion is being released at the inter- 
face of the solidified layer and the flowing warm liquid. The latent heat, along with the 
convective heat being transferred from the warm liquid boundary layer to the interface, 
is conducted through the frozen layer and wall and is then transferred to the coolant. 

The coolant also removes the additional amount of heat necessary to subcool the frozen 
layer and the wall as their temperatures decrease. The solid layer continues to grow 
until it achieves a steady-state thickness X_. Constant properties are assumed through- • 

out the analysis. 

I 


Steady-State Thickness of Frozen Layer 

In many solidification problems the frozen layer never approaches a steady-state 
thickness; rather, it continues to grow indefinitely with time. This occurs in a lake or 
river, for example, when the entire body of liquid is already at the freezing temperature 
and is exposed to a heat sink, the cold atmosphere in this instance. The only heat 
source in such a case is latent heat of fusion which is continually extracted so long as 
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the heat sink is present, thereby causing the frozen layer to grow continuously. 

When, however, the liquid flowing over the frozen layer is at a temperature above 
the freezing point, the warm liquid continually supplies heat to the frozen layer - liquid 
interface by convection. Under this condition the frozen layer eventually achieves a 
steady-state thickness. This thickness, which will be used as a reference length in the 
later analysis of growing layers, is derived here. 

It is assumed that the liquid-solid interface at the boundary of the frozen layer is 
always at the freezing temperature Tj. If the heat flow is taken as positive in the posi- 
tive x direction, the convected flux from the liquid to the liquid - frozen layer interface 
is at all times 


-q = hj (Tj - Tj) = constant 


At steady state a heat balance at the liquid - frozen layer interface gives the relation 


- q s - V T i 



— + — + — 
k kw h c 


This is solved for the steady-state thickness 


X 


s 




( 1 ) 


For X_ = 0, equation (1) gives the relation between variables required to just avoid 
s 

freezing. For example, solving for T^ gives 


x s =o 


= T f + 



( 2 ) 


To prevent freezing, the liquid temperature T^ must be equal to or greater than the 
value given by equation (2). 
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Description of Method of Analysis for Frozen Layer 
Growth and Temperature Distributions 


Before proceeding with the details of the analysis, it is instructive to outline the 
general features of the method. The heat flows within the frozen layer and the wall are 
governed by the one-dimensional transient heat conduction equations: 


k 




(3a) 


for the frozen layer and 


„ 8T W 

** ax 2 ' PwC I>w ar 


(3b) 


for the wall. Equations (3a) and (3b) are each integrated twice with respect to the space 
coordinate x. The first integration in each medium is between one of the boundaries and 
an arbitrary position x. This results in expressions for the local heat fluxes k(9T/9x) 
and k w (9T w /9x) in terms of the heat fluxes at the boundaries which are generally un- 
known. The second integration, again between one of the boundaries (not necessarily the 
same boundary used in the first integration) and an arbitrary position x, yields the local 
temperatures T(x, r) and T w (x, r) in terms of the temperature at one boundary of each 
medium and the heat flux at one boundary. 

In this problem the only known boundary temperature is at the liquid - frozen layer 
interface. The instantaneous temperatures and heat fluxes are unknown at the bounda- 
ries between the solidified layer and wall and between the wall and coolant. When the 
unknown boundary heat fluxes and temperatures are eliminated by some rather lengthy 
algebraic manipulations, there results two rather complicated coupled integral equa- 
tions, one each for T(x, t) and T w (x, r). 

When the equation for T(x, r) is evaluated at the liquid - frozen layer interface, 
which is at the known solidification temperature, an equation is obtained for the rate of 
frozen layer growth. This equation is integrated to obtain the layer growth, but the in- 
tegrated form contains integrals of T and T . The frozen layer growth equation and 
the two equations for T and T w form a set of three coupled integral equations. They 
are solved by an iterative method leading to closed-form approximate analytical solu- 
tions. 
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Temperature Distribution Equations 

i 

Equation (3a) is integrated from any position within the frozen layer x to the solid- 
liquid interface X: 


kil 

3x 


k^T 

9x 


= P c r 


,X 

/ ZE 
L 3t 


dx 


0 <x <X 


(4) 


At the liquid - frozen layer interface the heat conducted into the solidified layer is equal 
to that supplied by the latent heat of fusion and the convection from the flowing liquid: 


kil 

0X 


X 


= PL — + h (T. - T f ) 
dr L L 1 


(5) 


Equation (5) is substituted into equation (4) to give 


r x 

= +pL Z* + h, (T, - T f ) - pc / ZZ dx (6) 

dr 1 L 1 P 7 X 0r 

The term on the left side is the heat flow in the negative x direction crossing any posi- 
tion x at any time r. The last term on the right is the heat removed to subcool the 
portion of the solidified layer between x and X. The integration of equation (6) from 
the wall x = 0 to any position x results in an expression for the instantaneous temper- 
ature distribution in the frozen layer: 


k il 

3x 



where Tg = T 2 (t). 

In a similar fashion equation (3b) is integrated to obtain the temperature distribution 
in the wall. Integrating from the boundary in contact with the coolant x = -a to any x 
location within the wall gives 
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-a < x ^ 0 


( 8 ) 


V 


3T, 


w 

3x 


K 


3T. 


w 


3x 


= p c 

p. 


l-a 


w 


f X ^w 

L 


dx 


At x = -a the boundary condition is 


V 


0T. 


w 


SX 


- h c< T l - T c> 


-a 


where Tj = Tj(t). Substituting equation (9) into (8) gives 


V 


ST. 


w 

Sx 


" M T 1 " T c) P W C p w J 


ST. 


w 

St 


dx 


( 9 ) 

t 

i. 


( 10 ) 


To obtain the temperature distribution in the wall, equation (10) is integrated again from 
x = -a to x (the use of x = 0 to x as an alternate choice for the integration limits re- 
sults ultimately in the same final expression for T w (x, t)): 


T w (x, r) - Tj = — (Tj - T )(x + a) + 
w ^ 



-a < x s 0 


( 11 ) 


Equations (7) and (11) provide the temperature distributions in the frozen layer and 
wall but are not yet in usable form because they contain T^ and T 2 which are unknown 
functions of time. These temperatures can be eliminated by applying the conditions of 
continuity of temperature and heat flux that must hold at the boundary between the wall 
and frozen layer. Continuity of heat flux provides the relation 

k 9 ! 

3x 


x=0 


= V 


ST. 


w 


Sx 


x=0 


and by use of equations (6) and (10) evaluated at x = 0 this becomes 
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I 


r x r° 

PL f + h ! < T * - T f) - PC P J h c (T l - T C> + p w c p w J 

0 —21 


3T. 


w 

dr 


dx 


(12) 

Continuity of temperature requires that at x = 0, T w from equation (11) must equal Tgj 
therefore, 


T 2 - T i = — (T, - T )a + 
2 1 


HU 


0T. 


w 


0T 


dx Idx 


(13) 


Equation (12) yields an expression for Tj and this is substituted into equation (13) to 
obtain a relation for T 2 . The and Tg are then substituted into equations (7) and 
(11) to yield the temperature distributions 


T(x, t) = T c + ( 1 + 


V 


+ 1 (T f 

h„ dr h * 1 h„ A 


0T. 


c -a 


pc p f* 0T 
h c 4) 9r 


dx 


pL dX n l 


w 

3 T 


P c 


dx 


k dr k 


x + — (Tj - T f )x + 


w p. 


w 


- [U 


0T, 


^dxldx 

0T 


pc 

k 


■ftf 


3T 

0T 


dx Idx 


(14a) 


T w (x,r) = T c+ |l + — (x + a) | 


0T. 


PL dX + ^(T, f° 

h c dT h c h c /a 9t 


w 


dx 


„„ /.X 

^p f 

h A 


0T 

3t 


dx 


p c 
H w p 


w 




3X 


w 

3 T 


dxjdx (14b) 


L 



The temperature distributions in equation (14) can then be placed in dimensionless 
form by letting x> = x/X g , r' = h^T^ - T f )r/pLX g , X’ = X/X g , T’ = (T - T f )/(T C - T f ), 
= (T w - T f )/(T C - T f ), R = (X g /k)/[(l/h c ) + (a/k^)], S = c p (T f - T Q )/h and using 
equation (1) in the form 


x s h i < T i-T f ) _ R 
k (T f - T c ) 1 + R 

This gives 



A few transformations on equations (15) are now made. The time derivatives are taken 
out of the integrals by using the rule for differentiating under an integral. The double 
integrals are changed into single integrals; the method for accomplishing this is 
outlined in reference 1. Then the substitution (3/3 r’) = (dX’/dr' ) j( 3/ 3X* )J is made. This 
yields the temperature distributions in the forms 


T’(x’,X’) = 


/l + Rx’\/dX’ + 1 

\ RS dX' 

3I(x’,X’) ( n 

ay*)' 

V i + R / W’ , 

/ 1 + R dr’ 

3X ’ «w 

3X* 


(16a) 
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'1 + — Rx’> 


T^(x*,X’)= 1 - [ ^ )(^-+ l) -A 

\ 1 + R Adr’ ) 1 + 


where 


RS_ dX^ 
R dr’ 


aj(x',X') | a 9 J w ( x ' » X ' ) 

ax’ a... ax’ 


W 


(16b) 


i f x r x r* 

I( x * > X ) = i / T’ dx' + x* / T' dx' + / x'T' dx' 
R * / 0 Jx' JO 

VX ' > = Ti/a° T w® + /N T ;« 


J(x' 


f x ' rx' 

, x') = - / T' dx' + A X ' / T’ 
r y 0 kw J 0 


dx' 


(16c) 


(16d) 


(16e) 


= I T /a° T - + C * + * C <“ 


(16f) 


In the temperature distribution equations (eqs. (16a) and (16b)) there still remains 
the unknown quantity dX'/dr'. This quantity is also needed so that it can be integrated 
to obtain the frozen layer growth with time. The expression for dX'/dr’ can be found 
by imposing the physical boundary condition that at the liquid - frozen layer interface 
T = T f at x = X. This gives in equation (16a) T’ = 0 at x' = X', and equation (16a) can 
then be rearranged to yield 


where 


dX* 

dr’ 


R(1 - X') 


1 + RX’ + Rs(A_ + JL 

\dX' « w dX’ 


! rX' r X' 

G(X’ ) = I(x' = X’ , X' ) = — / T’ dx' + / x' T’ dx’ 

R-'O Jo 


(17a) 


(17b) 
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Equation (17a) is substituted into equations (16a) and (16b) to eliminate dX' /dr’ and 
to obtain the final forms of the temperature distributions: 


T» = R ( 1 ~ x> ) 
1 + R 


l + Rx’ + Rsf— + — 

\ 9X ’ « w dX«y 


1 + R 


R(1 - X’) 


1 + RX« + RS 


/ dG + 


\dX’ 


a dX' 
w 


(18a) 


R { 1 - — x’ 

V *w ) 


TL = 

W 1 + R 


1+ilW+Rsf— + — — ^ 

kw \ 3X ' 3X V 


1 + R 


R(1 - X’) 


1 + RX’ + RS ( — + — — ^ 
\dX’ a w dX’y 


(18b) 


Frozen Layer Growth 

Now return to equation (17a) and integrate to obtain the time variation of the frozen 
layer thickness. Separating the variables and integrating with the initial condition 
X’ = 0 at t’ = 0 yields 


X 1 X T X 1 

= f 1 + ft 3 ?... dX’ + S f — - — — dX' + S — f 
J 0 R(1 - X’) J 0 1 - X’ dX’ a w J 0 




dX’ 


1 - X’ dX’ 


(19) 


Integrating the first term of equation (19) directly and the last two terms by parts gives 


% 

i) 
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I 


< 


1 

♦ 



Equation (20) is the general expression that relates dimensionless frozen layer 
thickness and time. The quantities G(X') and I (X’ ) contain the temperature distribu- 
tions T'(x’ ,X') and T! (x’ ,X’) in the frozen layer and wall, respectively; these distribu- 
tions are provided by equations (18a) and (18b). The method of solution for the coupled 
temperature and layer growth relations will be given in the next section. The term 
1^(0) in equation (20) accounts for the temperature distribution in the wall at the instant 
that solidification begins. It should be noted that the influence of the wall heat capacity 
appears as the additive terms 



and 



in equation (20). 


Approximate Solution by Analytical Iterations 

The temperature distributions in equation (18) depend on the instantaneous frozen 
layer thickness, and the layer growth time as given by equation (20) depends on integrals 
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of the temperature distributions. An approximate solution to this coupled set of equa- 
tions is obtained by an iterative technique. 

The first approximations for the frozen layer and wall temperature distributions T' 
and and the growth time r' (X' ) are found by neglecting the effect of heat capacity in 

both media. When c and c are set equal to zero, the parameters S and l/o! w 

P P w 

become zero. Then from equations (18a), (18b), and (20) the first approximations are 

Tt = - x ') (21a) 

1 1 + RX’ 


and 




1 + RX’ 


(21b) 


Tt = -X* -i±Jhn(l -X’) 
1 R 


( 22 ) 


Improved approximations for T’, T^, and r’(X’) can now be obtained by substitut- 
ing equations (21a) and (21b) into the integral quantities G, I, I w , J, J w in equations 
(18a), (18b), and (20). When the integrations are carried out, the final equations for the 
second approximations take the form 


1 , T3 


1 + RX’ £ + RS 



1 + R 


R(1 - X’) 


dG, 


1 + RX’ + RS 


* + JL 


dX* a 


w 



(23a) 
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dG I _ R fx' + x ,2 + RX ,3 \ 

^ (1 + RX’) 2 \ R 3 / 

= (a* — - Ra' 2 + JL a.'A 

dX ' (1 + RX' ) 2 V k 3 \/ / 


£ = — and 0<|<1 
X 


T V,n^’ x ’) 



1 + R 



Ra'/3 + RS 





1 + R 


R(1 - X’) 


1 + RX’ + RS 


fdG 


dl 


idX' 


K-9L 


W T 




where dGj/dX’ and dl w /dX' are given by equations (23c) and (23d) and 


0J. 


l - U k 


ax* \R 


k (2 + RX’) 

Nv / 2 (1 + RX’) 2 


0 J, 


W T 


R 


dX ’ (1 + RX’) 2 ' R k 


^.a-Mia' 3 - a ’ 2 /3 2 _ R _k_ a ,3^3 






y3 = 25 and -1 < j3 < 0 
a 


and finally 


(23 c) 


(23d) 


(24a) 


(24b) 


(24c) 
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, 3 + 3R + R 2 ln (i _ x>) + ln(l + RX’) 


3(1 + R)' 


3R(1 + R)' 


_k_ 


X’R 


ln 


1 - X' 


(1 + R)(l + RX*) ^ + r ^2 1 + RX'J 


( 25 ) 


The solutions given by equations (23), (24), and (25) were the highest order approxi- 
mations that were carried out for the present problem. It might be possible to obtain the 
third order approximations, but the results would be so complex that they would not be 
useful for practical problems. A numerical approach would probably be simpler than 
carrying out the third approximations. The solution equation (25), however, can be 
readily evaluated and reveals the influence of the several independent parameters. In 

reference 1 where c was neglected, it was demonstrated that the second approximate 

P w 

solutions differed by only a few percent from the third approximations. It seems rea- 
sonable to assume that the same behavior holds true for the solutions presented in this 
report. 


Graphs for Predicting Solidified Layer Growth 


Equation (25) for predicting frozen layer growth can be written as 


T » = rt + SF + S — /a’ — - Ra’ 2 + -V 3 — 1W 

“wV k 3 V 


(26a) 


where rj, F, and W are functions of X’ and R only: 

1 + R 


Tj = -X’ - 


R 


ln(l - X’) 


(26b) 


F = -< 


RX’ L, + 2 + R 1 + 3 + 3R + R 2 ln(1 _ X1) + ln(l + RX’) 1 


3(1 + RX') L R(1 + R) 


3(1 + R)' 


3R(1 + R Y 


(26c) 
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1 


(26d) 


W = 


X'R 

(1 + R)(l + RX*) 


(1 + R)“ 


In 


1 - X* 

1 + RX* 


When the specific heat of the wall is zero, l/« w = 0, the contribution arising from wall 
heat capacity will be zero; that is, the term containing W vanishes. Similarly, SF is 
the contribution from the heat capacity of the frozen layer. The term Tj is the solution 
for no heat capacity in either wall or frozen layer. 

For convenience in making quick desk calculations, the quantities Tj, F, and W 
are given in figures 2(a), (b), and (c) as functions of R for various X' values. For a 
given set of imposed conditions (T c , , h c , h^ , a, etc. ) the R value is first computed. 
Then the dimensionless time r* to form a particular dimensionless thickness X* can 
be obtained by reading the quantities from the graphs and combining them according to 
equation (26a). This gives a relation between X' and t* from which X = X(r) can be 
found in dimensional form by using the definitions for X’ and t* . 


DISCUSSION 

As discussed in conjunction with the application of equation (26a), the three terms 
on the right side of the equation represent, respectively, the dimensionless growth time 
t* when all heat capacities are neglected, the contribution of the frozen layer capacity, 
and the contribution of the wall heat capacity. There are several parameters involved 
such as S, R, a/a w , and a’. A general parametric study of the effect of these inde- 
pendent parameters does not seem worthwhile, especially since equation (26a) can be 
easily evaluated for any design application. 

Before giving some illustrative examples, an examination of the last term on the 
right in equation (26a) (wall heat capacity term) will provide some additional information. 
The parameter R can be written in an alternate form: 


R = 


*Sv 


where 


y = 


1 


4 * 


h c a 
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For large h c a and low k^, y approaches 1 as a lower limit. On the other hand, for 
large k^ and small h c a, y can become very large. With this range of y in mind, the 
factor in the last term of equation (26a) can be examined: 


a’ ^ - Ra’ 2 + R2at3 — 
k 3 *w 



(27) 


As y — 1, $ — and when y — °°, $—1. Hence, the <E> is confined within rather nar- 
row limits of magnitude and is positive. Now the quantity W is examined to determine 
its sign. The term (X’R)/(1 + R)(l + RX T ) is positive. The argument 
(1 - X’)/(l + RX’) < 1; therefore, ln(l - X’)/(l + RX') < 0, and W is positive for all X’ 
and R. It is then evident that the contribution to t 1 by the wall heat capacity is posi- 
tive; that is, the growth time r for a given thickness X is increased because of the 
extra heat extraction needed to cool the wall. Since the factor F is positive, the same 
effect is true for the frozen layer capacity. This is what would be expected from intui- 
tive reasoning. 

To demonstrate the type of results obtained from the analysis a few illustrative ex- 
amples are now given. To choose some reasonable values for the parameters, one of the 
experimental tests from reference 2 is used as a starting comparison. In those tests, 
warm water flowed over an Inconel plate that was cooled on the opposite side by flowing 
chilled alcohol. For the data shown in figure 3 the Inconel plate was 3/16 inch 
(0. 476 cm) thick and a steady-state ice layer about 0. 4 inch (1 cm) thick was formed. 

The data is a little above the analytical prediction, but the agreement is satisfactory for 
engineering purposes. Unfortunately the data does not provide a good check on all as- 
pects of the theory because the heat capacity terms are relatively small for the thick- 
nesses tested. Additional data would be desirable having larger frozen layer and plate 
thicknesses. 

With the experimental test in figure 3 as a starting point, the predicted influence 
on the ice growth of changing some of the conditions is shown in figure 4. Figure 5 il- 
lustrates the corresponding temperature variations in the wall and frozen layer. 

Figure 4(a) differs from figure 3 only by having the wall thickness increased from 
3/16 to 1 inch (0. 476 to 2. 54 cm). The resulting steady-state ice layer is relatively thin. 
The growth curves reveal that the ice capacity has a negligible effect while the wall ca- 
pacity significantly slows the frozen layer formation. In figure 4(b) the coolant temper- 
ature is decreased from -40. 5° to -400° F (232. 9° to 33. 2° K) with the remaining condi- 
tions kept the same as in figure 4(a). Lowering the coolant temperature causes the 
steady-state ice layer to become quite thick. As shown by the growth curves, the ice 
capacity now has a significant effect throughout the layer growth; the wall capacity is 
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only important during the early portion of the transient (the temperature distributions of 
fig. 5(b) show it is during this period that the wall undergoes most of its cooling). In fig- 
ure 4(c) the conditions remain the same as for figure 4(b) except that the water tempera- 
ture is increased from 53. 5° to 150° F (285. 1° to 338. 7° K). This causes the ice layer 
to be relatively thin and provide a minor heat capacity effect; the effect of the wall ca- 
pacity is large since the wall is substantially cooled by the low temperature coolant. 

The solid lines in figure 5 show the transient temperature distributions in the wall 
and ice as computed from equations (23) and (24) for the conditions in figure 4. The 
curves in the ice layer terminate at the various x values equal to the instantaneous ice 
thicknesses at the times shown. 

The dashed curves in figure 5 are solutions to the one-dimensional transient heat 
conduction equation for the wall before the frozen layer begins to form. The solutions 
were carried out for the case where the wall was subjected to the water flow on one side 
and was initially at a uniform temperature equal to the water temperature. Then at time 
t = 0 the convective cooling condition was applied to the other side of the wall. The 
dashed curve for the largest f is the temperature distribution at the instant when the 
surface of the plate in contact with the water reaches 32° F (273. 2° K), at which time 
freezing is assumed to begin. At this r the dashed curve does not agree precisely with 
the solid curve for t = 0 (as calculated from eq. (24)). This brings us to some com- 
ments about the initial condition imposed on the freezing process. 

An instant before a frozen layer forms, r = 0”, the wall surface at x = 0 is at the 
freezing temperature of the liquid and the heat flux into the surface is that supplied by 
convection: 


T w (x = 0, r = 0") = T f 
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An instant later, r = 0 + , the frozen layer has just started to form and the heat balance at 
the surface becomes 
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A comparison of these equations for k w (0T w /0x) | x _q at r = 0” and r = 0 + reveals 
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that the heat flux into the surface at x = 0 undergoes a step increase in heat flux by an 
amount pL(dX/dT) | x _q when freezing begins. The effect of this additional heat flux is 
to retard the cooling of the wall in the vicinity of x = 0. In the remainder of the wall 
this effect has not been felt at small t and the temperatures continue to decrease in the 
same manner as before the freezing started. As a result, the temperature profiles 
change in shape from those shown by the dashed curves to those shown by solid curves. 

The present solution being approximate, since only the second approximation was 
carried out, does not reveal in precise detail the transient adjustment that occurs at 
small t. As shown by equation (20) the initial wall temperature distribution enters the 
successive approximations in an integrated form by means of the I w (0) term. The rea- 
sonably good agreement of the dashed and solid curves when the surface at x = 0 
reaches the freezing temperature shows that the freezing solution given here corre- 
sponds well with the following initial condition at the beginning of the entire transient 
process. When t < 0 there is no coolant flowing and the wall and flowing warm liquid 
are all isothermal at the temperature of the liquid that is flowing over one side of the 
wall. Then at f = 0 the coolant is introduced on the other side of the wall. 


CONCLUDING REMARKS 

In this report a solidification analysis is presented for conditions that arise in some 
important engineering applications such as casting of metals in molds, heat exchangers 
using cryogenics, and freezing of ice sheets on rivers. The problem studied is the tran- 
sient solidification of a flowing warm liquid in contact with a plane wall that is suddenly 
cooled convectively from the opposite side. The purpose of the analysis was to derive a 
method whereby the instantaneous frozen layer thickness could be predicted at any time, 
and which would account for the heat capacities in both the frozen layer and wall. Al- 
though this is a complex problem, a closed-form analytical solution (eq. (26)) was found 
that can be easily evaluated for engineering use. Portions of the equation are presented 
in graphical form to further facilitate its use. 

To demonstrate the type of results that are obtainable from this equation some ex- 
amples are given for ice forming on a stainless-steel plate. The examples illustrate 
that under some conditions the heat capacities of the frozen layer and wall appreciably 
retard the frozen layer growth. 

Lewis Research Center, 

National Aeronautics and Space Administration, 

Cleveland, Ohio, October 18, 1967, 

129-01-11-06-22. 
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Figure 1. -One-dimensional model for transient solidification of 
a flowing liquid on a convectively cooled wall. 



(a) Function xj. 


Figure 2 . - Graphical presentation of functions of X’ and R 
contained in equation (26) for frozen layer growth. 
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Wall capacity factor, W Frozen la y er ***** faclor ' 



(b) Relation between the frozen layer capacity factor F and the 
thermal resistance parameter R. 



Thermal resistance parameter, R ■ (X s /k)^[{l/h c ) + (a/k^ 

(c) Relation between the wall heat capacity factor W and the ther 
mal resistance parameter R. 


Figure 2. - Concluded. 






Time from start of ice growth, t, sec 


Figure 3. - Comparison of analysis with experimental data for 
ice forming on an Inconel plate 3/16 inch (0.476 cm) thick. 
Water temperature, Tp s 53. 5° F (285. l a K); coolant temper- 
ature, T c - -40. 5° F (232. 9° K) ; heat-transfer coefficients, 

h| * 92 Btu/(hrKft 2 l(°F) (522 W/(m 2 K°K)|, and h c * 118 
Btu/(hr)(ft 2 )(°F) (670 W/(m 2 )(°K)). 
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Frozen layer thickness, 



(b> Water temperature, 53.5° F (285. 1° K) ; coolant temperature, 



10° 10 1 10 2 10 3 10 4 
Time from start of ice growth, t, sec 


(c) Water temperature, 150° F (338. 7° K); coolant temperature, 
-400° F (33. 2° K). 

Figure 4. - Effect of heat capacity in prediction of ice layer 
growth on a thick plate. Heat-transfer coefficients: 

U t - 92 Btu/(hr)(fK)(°F) (522 W/(m2)(°K)) ; 
h c - 118 Btu/(hr)(ft 2 )(°F) (670 W/(m 2 )(°K)). 
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(a) Water temperature, 53. 5° F (285. 1° K); coolant temperature, 
-40. 5° F (232. 9° K). R - 0. 913, S = 0. 232. 
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(b) Water temperature, 53. 5° F (285. 1° K); coolant temperature, 
-400° F (33. 2° K). R = 10. 40, S * 1. 38. 
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(c) Water temperature, 150° F (338. 7° Kl ; coolant temperature, 
-400° F (33. 2° K). R = 1.08, S = 1.38. 


Figure 5. - Transient temperature distributions in wall and frozen layer. 
Heat -transfer coefficients: hj* 92 Btu/(hr)(ft2)(°F) (522 W/(m )( K)); 
h c = 118 Btu/(hr)(ft2)(°F) (670 W/(m^)(°K)>. 
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